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(57) Abstract: Image processing system for generating 
a multidimensional adaptive oriented filter to be 
applied to the point intensities of a d-dimensional 
image, comprising analyzing means with means 
(5, fi) to estimate at each image point a probability 
measure (Fi) of the presence of a type of feature of 
interest and a weighting control model (10) issuing 
a weighting control vector (11, VC) constructed 
from said probability measure, for the user to control 
synthesized adaptive kernels at each image point; and 
synthesizing means for generating the filter kernels at 
each image point adapted to the type of the features of 
interest, whose filtering strength is controlled by the 
weighting control vector. The system may comprise a 
selection unit (40) for the user to select synthesizing 
means for generating "pre -mixing filtering means" 
comprising combining means (30, XH) dependent 
on the type of the image features having inputs for 
the weighting control vector (11, VC) and the image 
data [I(x)] and having an output for weighted adaptive 
kernels (35, II) adapted to the type of the image 
features to produce the filtered image signal [H(x)]; 
and/ or "post-mixing filtering means" comprising 
both isotropic and anisotropic filtering means [15, gi)] 
applied independently of the type of the image features, 
whose outputs (Gi) are combined at each image point 
and adapted using the weighting control vector (11, 
VC) to produce the filtered image signal [G(x)]. 
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"IMAGE VIEWING SYSTEM AND METHOD FOR GENERATING FILTERS FOR 
FILTERING IMAGE FEATURES ACCORDING TO THEIR TYPE" 

Description 

Field of the Invention 

The invention relates to an image processing system and an image processing method 
to generate adaptive filters, whose supports are adapted to the type of features of interest in 
an image to be processed. The invention more particularly relates to such processing system 
and method that are applied to medical images. The invention also relates to a medical 
examination apparatus having such an image processing system and means for display. The 
invention finds an application in the industry of image processing systems and medical 
examination apparatuses. 

Background of the Invention 

A method to smooth and filter images is disclosed in the publication "A Framework 
for Anisotropic Adaptive Filtering and Analysis of Sequences and volume" in IEEE 
International Conference on Acoustics, Speech and Signal Processing (ICASP); San 
Francisco, March 23-26, 1992; p.469-472, vol.3, by H. Knutsson, L. Haglund and G.H. 
Granlund. The approach described in the cited publication makes use of quadrature and 
steerable filters having a set of different orientations including at least three different 
orientations in the case of 2-D image processing, and at least six different orientations in the 
case of 3-D image processing. The algorithm written from this approach requires the 
application of all the filters involved in the filtering of one image. The resulting algorithm 
complexity has been found to be high. So, the calculation time is particularly long. 

Summary of the Invention 

The present invention has for an object to provide an image viewing system having 
processing means to generate an adaptive filter for filtering an image comprising features of 
different type. By type of the image features, it is to be understood either anisotropic features 
such as features having elongated oriented shapes or isotropic features such as features rather 
showing circular or spherical shapes. In particular, it is an object of the invention to provide 
an image viewing system having analysis means for analyzing the content of a noisy image 
representing various objects and for determining the type of the features of interest of the 
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objects in the image; and having synthesis means for synthesizing filters whose supports are 
adapted to processing the different features according to their different types, for filtering 
noise or enhancing these features in the image. 

It is an object of the invention to provide analyzing means to estimate at each image 
5 point a probability measure of the presence of a feature of interest; then, from said probability 
measure, to determine a weighting control model for the user to control adaptive kernels at 
each image point. 

Such an image viewing system is claimed in claim 1. 

It is a further object of the invention to yield synthesizing means to provide the user 

1 0 with filtering means, among which: 

first filtering means called "pre-mixing filtering means" comprising combining means 
dependent on the type of the image features having inputs for the weighting control means 
and the image data and having an output for weighted adaptive kernels adapted to the type of 
the image features; and/or 

15 second filtering means called "post-mixing filtering means" comprising both 

isotropic and anisotropic filtering means applied independently of the type of the image 
features, whose outputs are combined at each image point and adapted using the weighting 
control means. 

A characteristic of the system of the invention is that at each point, the oriented filter 
20 means and the isotropic filter means are mixed according to the degree of feature type present 
at each point, such as for example the degree of anisotropic type and isotropic type of said 
point. Several problems are solved during those two phases of analysis and synthesis: 

The analysis is simplified however taking into account different kinds of 
information such as feature orientation and feature dimensions. The control means provided 
25 to the user allows him to fully specify the type of the features that he wishes to enhance. It 
allows him to fine-tune the filter in an extremely fast way, le. much faster than the actual 
filtering time. 

The solution to these problems is based on: 

TP 

The use of the control model based on the estimation of a vector (Fi,..F m ) of the 
30 pixel-level probabilities of the presence in the image of a number m of features of interest. 
This feature probability vector is captured by a set of m operators fi, . . ., f m . The design of 
those operators depends on the type of the features to be captured. 

The control model is completed by an explicit user-driven combination means X u of 
the feature probabilities into a k-dimension control vector (c l9 . . .,Ck) . This combination 
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specifies the user's wishes in terms of: "which feature or feature combination should be 
enhanced and with what strength". 

Selector means may be used to select filtering means chosen among two chains of 
proposed filtering means. The first filtering chain, called "pre-mixing filtering means", starts 
5 from the real-time image and produces an output signal Rh using a real-time adaptive- 
enhancing filter kernel H. The definition of this kernel H is specified by the above defined 
combined control vector (ci,. . .,Ck) using a combination operator Xh. The second filtering 
chain, called "post-mixing filtering means", comprises a combining operator X g , which 
combines the outputs of m predefined adaptive-enhancing filters g and produces an output 

10 signal Rg. Each of these m predefined adaptive-enhancing filters g is dedicated to the 
enhancement of a specific feature. Their combination is also driven by the control vector 
(ci,. . .,c k ) T via the combination operator X g . This chain allows extremely fast tuning of the 
filter by the user since the adaptive filters gi,. . ^g* are computed independently of the user 
tuning. Only the combination operator X g needs to be recomputed at every new control 

15 adjustment. 

The system of the invention has very powerful means to constitute adaptive filters for 
filtering features of different types such as anisotropic oriented features or isotropic areas. 
These filtering means show very low complexity of calculation, resulting in substantial 
computation saving. In addition, the synthesis means of the system of the invention gives to 

20 the user the choice between first filtering means and second filtering means whose 

complexities are different and whose computation loads are also different. Besides, the 
computations means according to the invention can be efficiently implemented using very 
simple techniques. Another advantage of the system of the invention is that the required 
computation time of the second filtering means is especially short. This system of the 

25 invention can constitute a very fast and general multidimensional adaptive filter. 

The second embodiment of the system of the invention can be implemented so that it 
is minimally demanding in terms of input/output, thus making a very good filter candidate for 
general purpose processor implementation. In fact, the system of the invention presents the 
advantage that it does not require any memory means. Since at the present time, processing 

30 means are very efficient for producing results of calculations, but are not very efficient for 
extracting data from memories, hence, the system of the invention, which only comprises 
calculation means, is very efficient. 
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In these advantageous embodiments of the invention, the user disposes of control 
means for interacting with the analysis means to select the weight to be applied to the 
filtering of the features of interest in the image and/or to the synthesis means. 

The invention may refer to image data in which each data point (pixel or voxel) is 
attributed a single scalar intensity value. This can be extended to multi-spectral data in which 
several values are attributed to each data point. For example, in colored images, one may 
refer to Red, Green and Blue components. Multidimensional data can be dealt with, such as 
one dimension, spatial 2-D, spatial 3-D, spatial multidimensional data plus intensity 
dimensional data and/or plus time dimensional data. When one of the dimension is time and 
when one aims at real-time results, the method of the invention may have recursive steps for 
using data recorded at previous instants. 

Brief Description of the Figures 

The invention is described hereafter in detail in reference to diagrammatic figures, 
wherein : 

FIG. 1 is a functional block diagram illustrating the means of the image viewing 

system; 

FIG.2 illustrates the steps carried out by the viewing system of the invention; 
FIG.3 is a functional block diagram of a viewing system and a medical apparatus 
according to the invention. 

Description of Embodiments 

The invention relates to an image viewing system having processing means to 
generate filters whose supports are adapted according to local type of features of interest in 
an image to be processed. The invention relates to such a system for adaptive filtering of 
multidimensional signals for quality improvement application including noise reduction and 
the enhancement of valuable details. This system has processing means for achieving this 
goal in two steps: first means for signal analysis and, using the results of this analysis, second 
means for filter synthesis. In particular, the invention relates to an image viewing system 
having analysis means for analyzing the content of a noisy image representing various 
objects, in order to determine the type of features of interest in an image; and having 
synthesis means for synthesizing filters with adapted supports apt to filter noise or enhance 
the features of different type in the image. The analysis means analyses the image signal in 
order to determine local properties of the image to be filtered. The synthesis means realizes a 
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filtering kernel, which varies at for each current point of the image and is adapted to the local 
properties of the image signal, while scanning the different points of the image during the 
filtering process. 

The system comprises means of acquisition of image data including point coordinates, 
5 denoted by x, and an intensity value, denoted by I(x), for each image point x. The proposed 
system deals with a multi-dimensional image, whose dimension is denoted by d. It is not 
excluded that d = 2. For a d-dimensional image, each data point is defined by a d-tuple point 
x associated to coordinates of the form (x l9 x 2 ,. . .x d ) where each of the components xi, 
x 2 ,. . .x d varies in a predefined range, for instance 0<xi<l! ; 0<x 2 <I 2 ; 0 <x d < I d . To each such 
10 d-tuple, corresponds an image value, which may be a single scalar value. In the case of a 
monochromatic image, the scalar value is the luminance intensity. For a color image, the 
image value may be represented by a set of three intensities representing the Red, Green and 
Blue components. 

FIG. 1 is a block diagram of the image processing system of the invention. Referring 
1 5 to FIG. 1 , the input image data I(x) are processed using the following processing means: 

1) Analysis means for analyzing the content of the original image. 
The Analysis means comprises a bank of a number m of operators, denoted by 
fi,..,fi,..,f m and referenced 5 in FIG.l. The bank 5 estimates at each current point of the image 
a probability of presence of features of interest. As feature of interest, one is to understand 
20 either an anisotropic feature which is an elongated oriented feature, such as for example in a 
displayed medical image representing vessels, an edge or a border of a vessel; or an isotropic 
feature which may be rather a circular or spherical feature, such as an aneurysm or a vessel 
perpendicular to the plane of the displayed image. 

More generally, for a 2-D image signal, anisotropic features might include crest, 
25 valleys, crest intersections, pass, plateau borders. It is to be noted that these typologies (crest, 
valleys, etc. . ..) refer to 2-D image signals where the intensity of the 2-D image is represented 
in a third dimension. Hence, the intensity is figured out as the height dimension in a 
landscape. Isotropic features might be "flat" isotropic areas, where flat means having constant 
intensity in the height dimension. For a 3-D image signal, the features might even be more 
30 diverse, including surfaces oriented in any directions, and 3-D planar or non-planar (skew) 
curves. 

The outputs of the bank of operators 5 are formed by a number m of probability 
values, denoted by Fi,..,Fi,..,F m , for the m types of different features to be present in the 
image to be filtered. 
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In an example, in an image, two types of features may be present, such as edges and 
isotropic areas. The probability of presence of edges is expressed by a, and the probability of 
presence of isotropic areas is expressed by (1 - a). 

In another example, in an image, three types of features may be present, such as 
5 edges, edge intersections, and isotropic areas. The probability of presence of simple edges is 
expressed by a (alpha), and the probability of presence of isotropic areas is expressed by p 
(beta), and the probability of presence of edge intersections is expressed by y (gamma). The 
total probability value is preferably normalized so that: 

a + p+y=l. (1) 

10 The general formula giving the probability of presence of different features may be expressed 
by: 

£FW (2) 

i=l 

The type of the operators of the bank 5 depends on the type of the features of interest. 
It is not the object of the invention to specify these types very strictly. In the above described 
15 case when edges and isotropic areas may be found as features in the image, the probability a 
of presence of edges can be determined for example by calculating the intensity gradient at 
each point in the image. 

Hence, according to the invention, a number m of features of different type are found 
in the image using a corresponding number m of operators 5 that depends on the type of the 
20 features. 

A combining operator 10, called user operator, denoted by X u , further receives at its 
input the m probabilities F, of presence of features of different type and provides at its 
outputs a control vector of k components Ci,..,Q,..,Ck. This control vector 11 of k 
components, denoted by V c , allows to controlling the adaptive kernels of the filter to be 

25 constructed according to the invention. The number k of vector components is such that m ^ 
k, depending on whether the user wants to group the filtering of different types of features 
together or not. Using this control vector V c , the user will be able to define the feature or the 
group of features that he wants to filter or to enhance and the strength of the filtering or 
enhancement. For instance, in the above-described example of an image comprising edges 

30 and circular features, the user may choose to enhance the isotropic features with strength of 
0.8 and the anisotropic features with strength of 0.2. 

Hence, the analyzing means provides the control vector V c , which is apt to control 
the adaptive kernel to be passed on each point in the image. It is to be noted that this control 
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vector depends on the type of the features and not on the direction of particular features such 
as the direction of anisotropic features. 

2) Synthesis means for synthesizing oriented filters. 

The synthesis means comprises two possible chains of means for generating an 
5 adaptive filter. The first chain is denoted by H and called "pre-mixing filtering chain of 
generating adaptive filter". The second chain is denoted by G and called "post-mixing 
filtering chain of generating adaptive filer". The user has means 158 called user control 
interface, for selecting the first or the second chain. 

2.1) Pre-mixing filtering chain of generating adaptive filter. 
10 This chain comprises a combining operator denoted by Xh, referred to as 30, which 

receives at its input, the control vector V c and the image data I(x) and which provides at its 
output an adaptive kernel. This adaptive kernel is adapted to the orientation of the anisotropic 
oriented features and/or to the dimensions of the isotropic features to be filtered or enhanced. 
In addition, this adaptive kernel is adapted to respect the weighting parameters defined by the 
15 control vector Vc, since Vc is used on top of orientations and/or dimensions of the features. 

Those skilled in the art know means to synthesize a filter defined by a probability 
such as a and by structure tensors. A map of the a probabilities defines the moments of a 
target tensor. From the target tensor, a filter is synthesized where the tensor eigen vectors 
give the local orientation and where the tensor eigen values give indications of the strength in 
20 the different orientations. Hence, the information of strength along different directions that 
are given by the control vector may be used to control a tensor. According to the invention, a 
tensor is synthesized, in order to detect anisotropic or isotropic features, so that it is aligned 
on a local tangential direction. 

In an example, the number m of features of the image is m = 2, the first feature 
25 correspond to Fi, and the second feature to F 2 . The tensor is defined by a tangential moment 
\i t and by a perpendicular moment (i p . The minimum and maximum moment are respectively 
Umin and Umax. The moment are given by the following formulae: 

M* = Umin + <* (M* max - M-min ) (3a) 
Hp = Umin + cc (jj, p max - Umin )• (3b) 
30 When a = 0, \i t = \i P = I^min ( 4 ) 

hence the tensor is isotropic since the strength is the same in the directions t and p. 

When a = 1, \it = |a t max, ( 5a ) 311(1 

M-p = Hp max ( 5b ) 
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hence the tensor is completely oriented. 

With the above-defined tensor, the filter can be realized: this filter will clearly depend 
on the type of the features that are present in the image and that the user wants to favor. In 
this chain of filter generation, the information about the type of features is mixed to perform 
5 the generation of the filter kernel H referred to as 35. The result 3 of the filtering is given by 
H(x). 

2.2) Post-mixing filtering chain of generating adaptive filter. 
In this second chain of filter generation, instead, the image data are first filtered using 
a bank 15 of a number m of different filtering means denoted by gi,.->giv->gm- The number of 
10 these filtering means corresponds to the number m of features of different type to be 

processed. The outputs, denoted by Gi,..,Gi,..,G m , of the different filtering means are mixed. 
Each different filtering operation relates to a given identified feature. In this chain, adaptive 
or non-adaptive filters are used. Each of these filters is specialized in the filtering of one 
given type of feature, that has been previously identified, such as edges, intersections and 
1 5 circular or spherical features. 

The combination of the outputs of the filters of the bank 15 is performed using a 
combination operator Xq. This combination operator Xg, referred to as 20, also receives the 
control vector Vc of the k components Ci to Ck. 

For instance, the combination may be a simple weighted sum of the results of the 
20 filter bank 15. This combination operator Xg is very important element of this chain, because 
it allows to changing very easily the weights related to the different type of features in the 
final result 2 denoted by G(x), whichever are the filtering means of the bank of filters 15. The 
filtering means of the bank 15 are applied independently of the types of the features 
examined in the image. The filtering operations using the bank 1 5 are realized in a 
25 preliminary step and the weight variations are performed while combining the results 
Gi,..,Gi v .,G m , by using the control vector Vc. 

In an example, an image comprising anisotropic features and isotropic features can be 
filtered using two types of filters denoted by gi and g 2 in the bank of filters 15. Gi is for 
example the output of the oriented filter gi and G 2 is the output of the isotropic filter g 2 . The 
30 filter g2 is not necessarily adaptive, it may be a simple Gaussian kernel. The user may choose 
a filtering strength for the anisotropic features given by a. In this case, the result at the output 
of the combining operator X G , referred to as 20, is given by the formula: 

G(x) = cc(x) Gi(x) + (1 - a) G 2 (x) (6) 
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Since the value a(x) depends on the choice of the user, the user can change the result 
G(x) according to its desire to improve the filtering strength or to improve the enhancement 
of the different types of features. This operation of varying the filtering strength can be 
achieved without performing any new calculation of the results Gi,..,Gj 9 ..,G m of the different 

5 filters gi,..,gt,..,gm of the bank of filters 15. 

In an example of realization, the oriented filter gi can be made using structure tensors 

with: 

jit = a constant (7a), and 

fi p = another constant. (7b) 
10 It is to be noted that the isotropic filter g 2 does not depend on (J* or m, values. 

3) Image processing method. 

Referring to FIG.2, the image processing method, for generating a multidimensional 
adaptive oriented filter to process image data in a number d of dimensions, operated by the 
above-described system comprises steps of: 

15 analyzing (5, £) the image to estimate at each image point a probability measure (Fj) 

of the presence of a type of feature of interest and to determine from said probability measure 
a weighting control model (10) yielding a weighting control vector (11, Vc) for the user to 
control synthesized adaptive kernels at each image point; and 

synthesizing filter kernels at each image point adapted to the type of the features of 

20 interest, which filter kernels are controlled by the weighting control vector. 

The user may select a chain of steps for generating the filter kernels called "pre- 
mixing filtering chain of steps", comprising a combination operation (30, Xh) dependent on 
the type of the image features for combining the weighting control vector (1 1, Vc) and the 
image data [I(x)] and yielding weighted adaptive kernels (35, H) adapted to the type of the 

25 image features to produce the filtered image signal [H(x)]. 

The user may select another chain of steps for generating the filter kernels called 
"post-mixing filtering chain of steps", comprising performing different types of filtering 
operations [15, ©)], such as isotropic and anisotropic filtering operations, applied 
independently of the type of the image features. These filtering operations yield signals (GO 

30 that are further combined, at each image point of the image, using the weighting control 

vector (11, V c ), in order to produce the filtered image signal [G(x)] adapted to the types of 
features. 

4) Medical viewing system and apparatus 
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The above-described means and steps can be implemented in the viewing system of 
the invention. Fig. 3 shows the basic components of an embodiment of an image viewing 
system in accordance to the present invention, incorporated in a medical examination 
apparatus. The medical examination apparatus 151 may include a bed 1 10 on which the 

5 patient lies or another element for localizing the patient relative to the imaging apparatus. 

The medical imaging apparatus 151 may be a CT scanner or other medical imaging apparatus 
such as x-rays or ultrasound apparatus. The image data produced by the apparatus 151 is fed 
to data processing means 153, such as a general-purpose computer, that carries out the steps 
of the method. The data processing means 153 is typically associated with a visualization 

10 device, such as a monitor 154, and an input device 155, such as a keyboard, or a mouse 156, 
pointing device, etc. operative by the user so that he can interact with the system. The data 
processing device 153 is programmed to implement a method of processing medical image 
data according to invention. In particular, the data processing device 153 has computing 
means and memory means to perform the steps of the method. A computer program product 

15 having pre-programmed instructions to carry out the method may also be implemented. 

The drawings and their description herein before illustrate rather than limit the 
invention. It will be evident that there are numerous alternatives that fall within the scope of 
the appended claims. Moreover, although the present invention has been described in terms 
of generating image data for display, the present invention is intended to cover substantially 

20 any form of visualization of the image data including, but not limited to, display on a display 
device, and printing. Any reference sign in a claim should not be construed as limiting the 
claim. 
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Claims 

1 . Image processing system for generating a multidimensional adaptive oriented filter to 
be applied to the point intensities of an image formed in a number d of dimensions, 
comprising: 

5 analyzing means comprising means (5, fi) to estimate at each image point a 

probability measure (Fi) of the presence of a type of feature of interest and to determine from 
said probability measure a weighting control model (10) yielding a weighting control vector 
(11, Vc) for the user to control synthesized adaptive kernels at each image point; and 

synthesizing means for generating filter kernels at each image point adapted to the 
10 type of the features of interest, which filter kernels are controlled by the weighting control 
vector. 

2. The image processing system of Claim 1 , wherein the synthesizing means comprises 
means for generating: 

15 filtering means called "pre-mixing filtering means" comprising combining means (30, 

Xh) dependent on the type of the image features having inputs for the weighting control 
vector (1 1 , V c ) and the image data [I(x)] and having an output for weighted adaptive kernels 
(35, H) adapted to the type of the image features to produce the filtered image signal [H(x)]. 

20 3. The image processing system of one of Claims 1 or 2, wherein the synthesizing 
means comprises means for generating: 

filtering means called "post-mixing filtering means" comprising both isotropic and 
anisotropic filtering means [15, gi)] applied independently of the type of the image features, 
whose outputs (Gi) are combined at each image point and adapted using the weighting 

25 control vector (11, Vc) to produce the filtered image signal [G(x)]. 

4. The image processing system of one of Claims 1 to 3, wherein the analyzing means 
comprises a number m of operators (fi,..,fi,..,f m ), which outputs at each current point of the 
image a probability value (Fi,..,Fi,..,F m ) of presence of features of interest among the m types 
30 of different features in the image to be filtered. 



5. The image processing system of one of Claims 1 to 4, wherein the analyzing means 
comprises a combining operator (10, X u ), called user operator, which receives at its input the 
m probabilities (F lv .,F iv .,F m ) of presence of different types of features for forming the 
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weighting control model, which provides at its outputs the control vector (V c ) of k 
components (Ci,..,Ci,..,Ck), for controlling the adaptive kernels of the synthesized filters. 

6. The image processing system of one of Claims 1 to 5, wherein the features of interest 
5 are chosen among anisotropic features and isotropic features. 

7. The image processing system of Claim 6, wherein the synthesizing means comprises, 
in the "pre-mixing filtering means", a combining operator (30, X H ), which receives at its 
input, the control vector V c and the image data I(x) and which provides at its output an 

10 adaptive kernel (H) that is adapted to the orientation of the anisotropic oriented features 

and/or to the dimensions of the isotropic features to be filtered or enhanced by the weighting 
parameters defined by the control vector (V c ). 

8. The image processing system of Claim 6, wherein the synthesizing means comprises, 
15 in the "post-mixing filtering means", a number m of different filtering means (gi,..,gi,.->gm)» 

which are applied independently of the types of the features examined in the image, whose 
number m corresponds to the number m of features of different types to be processed and 
whose outputs (Gi,..,Gi,..,G m ) are mixed in a combination operator (Xg), which is controlled 
by the control vector (V c) to produce the filtered image signal [G(x)] . 

20 

9. The image processing system of Claim 8, wherein the combination operator (X G ) is a 
weighted sum of the results of the different filtering means (gi,.>>gi>">gm)» 

10. The image processing system of one of Claims 2 to 9, wherein a user control interface 
25 (158) is provided for the user to control the weighting parameters (Ci,..,Ci,..,C k ) of the 

control vector (V c ); for selecting the "pre-mixing filtering means" and/or the "post-mixing 
filtering means"; for varying the strength of filtering or enhancement related to the different 
type of features at the output of the combination operator (X G ) while combining the results 
(Gi,..,Gi,..,G m ) to produce the filtered image signal [G(x)] in the "post-mixing filtering 
30 means"; and/or at the input of the combination operator (X H ) in the "pre-mixing filtering 
means". 

1 1 . Image processing method for generating a multidimensional adaptive oriented filter to 
process image data, using a system as Claimed in one of Claims 1 to 10. 
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12. Medical examination apparatus comprising means to acquire d-dimensional image 
data [I(x)], a system as Claimed in one of Claims 1 to 10 and further comprising a display 
system (154) for visualizing processed images and user control means (158) for selecting 

5 weighting parameters and/ or acting on the user operator (10) and/or the selection unit (40). 

13. A computer program product comprising a set of instructions for carrying out the 
method as claimed in Claim 1 1 . 
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